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£N| ■ Abstract. We present here a complete list of contributions to the gravita- 

tional wave frequency evolution from compact binaries on circular orbits up to 
the second post-Newtonian order by including the interaction of magnetic dipole 
■ moments, quadrupole-monopole interactions together with the spin-orbit, the 

spin-spin and the self-interaction spin contributions. We apply our results to 
the Manchester et al. model of the J0737-3039A-B double pulsar. 
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\0 '. 1- Introduction 

o 

Inspiralling compact binaries (like double-pulsar systems) emit gravitational 
radiation. These binary systems are the most promising source of the gravi- 
tational waves for the Earth-based interferometric detectors (LIGO, VIRGO, 
GEO, TAMA and AIGO). The coalescence of such compact binaries is preceded 
by a milder inspiral phase, for which the post-Newtonian approach provides a 
reliable description. This description is generally considered valid until the sys- 
tem reaches the innermost stable circular orbit (ISCO). The ISCO is not even 
precisely defined. For the special case of a non-rotating black hole the ISCO is 
r=6M and at r=M for a maximally rotating black hole (throughout the paper we 
use G = c = 1). For neut ron stars the gravitational wave frequency is 800-1230 
Hz (|Oechslin et al. 2004ft at ISCO. The LIGO sensitivity is 1000 Hz. 
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The description of motion was given to 3.5 post -Newtonian (3.5 PN) order 

accuracy, with the in clusion of spin- orbit effects fSO \dWex 19951: iBieth fc Schafer 19971 : 
Ger gelv et al. 1998) an d their first P N correction (Ta goshi et al. 200lD. Spin- 
spin fSSl dKidder 199.51: iGergelv 200 01 auadrupole-monopole (OM) dPoisson 19 98 : 
Gergelv fc Keresztes 2003ft~ and magnetic dipole-magnetic d ipole (DDI (lloka <fe Taniguchi 2000 ; 



Vas uth et al. 2003ft and"self-interaction spin term (SS-self) (|Gergelv 20 00: Mi koczi et al. 2 005) 



contributions to the accelerations were also discussed. On the long run due to 
the emission of gravitational waves the orbit tends to circularize. Therefore we 
consider circular orbits, for which the gravitational wave frequency is twice the 
orbital frequency. We evaluate the rate of increase of f. This is given by the 
rate of change of the orbital angular frequency uj = irf under radiation reaction. 
We present the accumulated number of cycles left until the ISCO and the spin 
parameters in the Manchester et al. model of the J0737-3039A-B double pulsar. 
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2. Frequency Evolution and Accumulated Number of Cycles 

In this section we summarize some of the results of (|Mikoczi et al. 2 005) on 
the frequency evolution and accumulated number of cycles for compact binaries 
on circular orbit. The radial projection of the acceleration defines the orbital 
angular velocity ui. We expressed both the energy and the energy loss with oj. 
Then the evolution of the radiative orbital angular frequency is: 
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0" = 4 ^SS—self 4 &QM 4 &DD- (2) 

The quantities /3, crs 1 s 2 , <?SS-self> a QM an d o"dd are the spin-orbit, spin-spin, 
self-interaction spin, quadrupole-monopole and magnetic dipole-dipole parame- 
ters, respectively: 
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The N, P N, SO, SS, 2PN and ta il contributions were verified to a gree with thos e 



given in (P oisson W ill 199 51), the QM with tho se given in ((Poisson 1 998). 



the DD with those given in flloka fc Tanitruchi 20 00) and the SS-self with those 
given in (iMikoczi et al. 20051 ). respectively. The constant Aq is also given in 



( Vasuth et al. 20031 ). 



From here the accumulated number of gravitational wave cycles emerges as: 
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where the dimensionless time variable r = rj(t c — t)/5m is related to the time 
(t c — t) left until th e final coalescence. T he tail and 2PN contributions agree 
with those given in dBlanchet et al. 2002ft . 



Tabic 1. Accumulated number of gravitational wave cycles. 



PN Order 


J0737 - 3039 


51913 + 16 


BH -BE 




1.337M Q 


1.387M© 


10 4 M Q 




1.250M Q 


1.441M© 


10 5 M Q 


fin(Hz) 


10 


10 


4.199 x 10" 4 


ff in (Hz) 


1000 


1000 


3.997 x 10~ 2 


N 


18310 


15772.1 


21058 


PN 


475.8 


435 


677 


SO 


17.5(3 


16.5/3 


36/3 


SS - self, SS, QM, DD 


-2.1(7 


-2.1(7 


-5(7 


Tail 


-208 


-206 


-450 


2PN 


9.8 


9.5 


18 



Table 2. The spin parameters for the two solutions (JR1 & JR2) of the 
Jenet-Ransom model and for two solutions (Ml and M2) of the model of the 
Manchester et al. representing the binary pulsar J0737-3039. 



Spin parameters 


JR1 


JR2 


Ml 


M2 




[k x = 167°] 


[ki = 90°] 


[«i = 14°] 


[ki = 60°] 




-0.166 


0.001 


0.168 


0.087 


°s lS2 (10- 4 ) 


-0.372 





0.371 


0.191 


(?SS-self (10 4 ) 


0.298 


0.345 


0.298 


0.333 



In Table 1 we present numerical resul ts for two well-known binary neutron 
star systems, the double pulsar J0737-3039 ( Burgav et al. 200*31 : iLvne et al. 20041 ) . 
Hulse- Taylor pulsar B1913+16 and one example of galactic black hole-galactic 
black hole binary. 



3. Application for the J0737-3039A-B Double Pulsar 

The double-pulsar J0737-3039 has become an important astrophysical laboratory 
for testing gravitational physics. The measurements on this system where done 
by th e Parkes 64-m radio-tel escope at three frequencies: 680, 1390 and 3030 
MHz (|Manchester et al. 2004ft . 

We evaluate the spin parameters (Ts 1 s 2 an d &ss-seif for the case of the dou- 
ble pulsar J0737-3039. The neutron stars in this double pulsar have average radii 
15 km, pulse periods of 22.7 ms and 2773.5 ms. Tha k\ a ngle is between the spin 
axis and angular momentum. The Jenet-Ransom model ( Jenet k, Ransom 2004ft 



has two possible values: K\ = 167° ± 10° (JR1), and k\ = 90° ± 10° for pulsar A 
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(JR2). The other angle K2 for pulsar B in principle can b e determined by so lving 
numerically the spin precession equations. According to ( Kaspi et al. 2004h it is 
likely that wind-torques from the energetically dominant component have driven 
the spin axis of the other component to align with the direction of L, causing 
K2 = 0. Therefore 7 = k\. We give the estimates for the spin parameters in 
Table 2. 

Besid es the JR1 and the JR2 solutions we present here for the first time 
an model ( Manchester et al. 20041 1. According to this model the measurements 
the favorize maximally allowed the angle k\ = 60°. The best-fit solution (Ml) 
is located at m = 14°. We also present results for the maximally allowed value 
of K\ (M2). As it can be seen in this model the proper spin-spin contribution is 
comparable with the self-interaction spin contribution. 



4. Conclusions 

We have presented the complete set of contributions up to second post-Newtonian 
order (PN, SO, SS, QM, DD, tail and 2PN) to the evolution of gravitational wave 
frequency and to the accumulated number of gravitational wave cycles left until 
the final coalescence, with the inclusion the self- interaction spin terms (SS-self). 

We have shown that on the Manchester et al. model of the J0737-30309A-B 
double pulsar the self-interaction spin contributions are comparable with spin- 
spin contributions. 

This work was supported by OTKA no. TS044665 and T046939 grants. 
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